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ABSTRACT 

CoRoT, the first space-based transit search, provides ultra-high precision light curves with continuous time-sampling over periods of up to 5 
months. This allows the detection of transiting planets with relatively long periods, and the simultaneous study of the host star's photometric 
variability. In this letter, we report on the discovery of the transiting giant planet CoRoT-Exo-4b and use the CoRoT light curve to perform a 
detailed analysis of the transit and to determine the stellar rotation period. The CoRoT light curve was pre-processed to remove outliers and 
correct for orbital residuals and artefacts due to hot pixels on the detector After removing stellar variability around each transit, the transit light 
curve was analysed to determine the transit parameters. A discrete auto-correlation function method was used to derive the rotation period of the 
star from the out-of-transit light curve. We derive periods for the planet's orbit and star's rotation of 9.20205 ± 0.00037 and 8.87 ± 1.12 days 
respectively, consistent with a synchronised system. We also derive the inclination, ; = 90.00^[! !!!!j in degrees, the ratio of the orbital distance to the 



stellar radius, a/R^ = 17.36^jJ"25, and the planet to star radius ratio Rp/Rg 



l047+o'"'« 



We discuss briefly the coincidence between the orbital 



period of the planet and the stellar rotation period and its possible implications for the system's migration and star-planet interaction history. 
Key words, planetary systems - techniques: photometry 



1. Introduction 

Transits provide unique insights into fundamental aspects of 
extra-solar planets that are currently beyond the reach of other 
techniques: mean density (through the determination of true 
masses and radii), atmospheres (through transmission spec- 
troscopy and secondary transit observations), and formation and 
evolution mechanisms (through the statistics of the orbital pa- 
rameters, including the true mass). All but two of the transit- 
ing planets published so far have been very short-period planets 
(< 5 d), whose properties are bound to be affected by extreme 
proximity to their host star. The except ions are HP 14 7506b 
(iBakos et al.| [2007) and HD 17156b (Barb ieriet al]|2007l) . with 
periods of 5.6 and 21. 2d, both in very eccentric orbits (e > 0.5). 

CoRoT is the first space-based transit survey, and over its 
lifetime will survey 120000 stars for up to 5 months with pre- 
cisions down to 0. 1 mmag per hour. This letter is the 4* in a 
series presenting new extra-solar planets discovered thanks to 
the CoRoT observ ations and associat ed grou nd-based foUow-u 
program. Papers I jBarge et al.ll2008l) and II (lAlonso et al.ll200 
present the discoveries of CoRoT-Exo-lb and CoRoT-Exo-2b , 
both very close-in giant planets. Paper III (iBouchv et al.ll2008l) 
reports on the observation of the spectroscopic transit of CoRoT- 
Exo-2b. Here we report on the detection of CoRoT-Exo-4b, a 



* The CoRoT space mission, launched on December 27* 2006, has 
been developed and is operated by CNES, with the contribution of 
Austria, Belgium, Brasil, ESA, Germany and Spain. The first CoRoT 
data will be available to the public from February 2009 through the 
CoRoT archive: http://idoc-corot ■ ias.u-psud. fr/. 



relatively long period (9.2 d) gas-giant transiting planet, and on 
the analysis of the CoRoT light curve (the production of which 
is described in Section |2]i to determine the transit parameters 
(see Section [3]l and the stellar rotation period (see Section |4]i. 
This analysis makes no assumptions aside from the planetary 
nature of the companion, which was confirmed through ground- 
based follow-up observations which are reported in a compan- 
ion letter (Moutou et al. 2008 hereafter Paper V). In Paper V, the 
full system parameters are determined using the transit param- 
eters reported in the present work. These parameters are used 
in Section |5] where we discuss the implications of the rela- 
tively long period and of its proximity to the stellar rotation pe- 
riod. A further dete ction, CoRoT-Exo-3b, will be presented in 
iDeleuil et alJ (l2008l Paper VI). 



2. Observations and data processing 

Like the previously published CoRoT planets, CoRoT-Exo-4b 
was first detected during near real-time analysis of raw data, the 
'alarm mode' (which is described in Papers I and II), trigger- 
ing the ground-based follow-up described in Paper V, which in- 
cluded ground-based photometry in- and out-of-transit, at higher 
spatial resolution than CoRoT's own, to identify which of the 
stars falling in the CoRoT aperture was being eclipsed, and 
multiple radial velocity measurements to derive the companion 
mass. Once the planetary nature of the companion was con- 
firmed, a high resolution, high signal-to-noise spectrum of the 
host star was obtained to derive accurate stellar parameters. 
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CoRoT-Exo-4 (GSC designation 0480002187), whose coor- 
dinates and magnitude are given in Table[Tl was observed as part 
of CoRoT's initial run, during which ~ 12 000 stars with magni- 
tude 12 < 7? < 16 falling in a 1.3° x 2.6° pointing close to the 
anticentre of the Galaxy were monitored nearly continuously for 
58 days, starting on the 6* of Februai-y 2007. A total of 72319 
flux measurements were obtained for CoRoT-Exo-4. For the first 
33 d of the run, the time sampling is 512 s, after which it was 
switched to 32 s as the transits were detected by the alarm mode. 

Aperture photometry through a mask, automatically selected 
from a set of 256 templates at the beginning of the run, is per- 
formed on board. For stars brighter than R = 15, the flux is split 
along detector column boundaries into broad-band red, green 
and blue channels. Although the transits were detected in the 
raw data, the analysis presented here was based on the pipe line - 
processed light curve. The pipeline dAuverpne et al.ll2b08h cur- 
rently includes background subtraction and partial jitter correc- 
tion. For each exposure, a global background level is estimated 
from a handful of 10 x 10 pixel background windows distributed 
over each CCD, excluding background windows affected by hot 
pixels, and subtracted. The jitter correction is based on the satel- 
lite line of sight information, which is derived from the aster- 
oseismology channel (which lies next to the exoplanet channel 
on the focal plane). Currently, the pipeline only applies a rela- 
tive jitter correction for the three colour channels. This correc- 
tion conserves the total (white) flux, and no jitter correction for 
the total flux is attempted. The pipeline also flags data points 
collected during the SAA or affected by other events likely to 
impair the data quality, such as entrance into and exit from the 
Earth's shadow. 

The blue channel light curve is affected by a hot pixel event 
a few days after the start of the observations, causing a sudden 
rise in the measured flux followed by a gradual decay of ~ 5 
days. This was corrected by fitting an exponential curve to the 
decaying segment, excluding a small portion that falls in-transit 
(see Figure[Tl online only). That section of the blue channel light 
curve also shows a gradual rise in flux. To preserve it, a linear 
rise between the local mean flux levels before and after the sec- 
tion of light curve affected by the hot pixel was added after sub- 
tracting the exponential decay. The red and green channels are 
free of visible hot pixel events. Note that the transit depth is the 
same in all three channels (within the uncertainties), as expected 
for a planetary transit. 

A single band-pass is sufficient for the present work, so 
the flux from the three colour channels was summed to give 
a 'white' light curve (approximately covering the range 300- 
1000 nm). The three-colour photometry will be discussed in 
an upcoming paper pending improvements in the pipeline. 
Additionally, a version of the light curve with regular 512 s time 
sampling was computed by rebinning there over-sampled part 
of the original. This version was used to study the out-of-transit 
variability, while the over-sampled version was retained to es- 
timate the transit paramete rs. A short-bas eline (5 data points) 
iterative non-linear filter ( Aigrain & Irwinli2004 ) with 5-cr clip- 
ping was applied to both the over-sampled and regularly sampled 
light curves to identify and reject further outliers, resulting in a 
final duty cycle of 87%. 

Based on prel iminary ground-based imaging of the field 
dPeleuil et alj|2006 ) and a basic empirical model of the CoRoT 
PSF, we estimate the contamination of the photometric aperture 
by other stars than CoRoT-Exo-4 to be 0.3 + 0.1% (the uncer- 
tainty arises from the PSF model). This estimate is compatible 
with the transit depths measured from the ground (see Paper V). 
We therefore subtracted a constant equal to 0.3% of the median 




Fig. 2. Pre-processed white light curve of CoRoT-Exo-4. The 
light curve was normalised by dividing it by its own median, 
which gives more weight to the later (oversampled) part of the 
light curve. This does not affect the transit analysis, since it is 
the local out-of-transit flux around each transit that is taken as a 
reference (see Section[3]l. 

Table 1. Star and transit parameters derived from the CoRoT 
light curve. The quantity mI^^ /R^, which is used in determin- 
ing the stellar parameters (see Paper V), is derived directly from 
a/Rs and P. 



Parameter 


Value 


Bayesian range 


RA 


06 48 46.70 




Dec 


-00 40 21.97 




R - mag 


13.45 




Prot (d) 


8.87 ± 1.12 





P (d) 9.20205 ± 0.00037 

To (HJD) 2454141.36416 ± 0.00089 

17.36+ 



90 000+° '""' 



a/R, 
u 

mI'^/r, 



r+O.OS 

-0.25 



44+0.16 

1047+°''°4i 
899+0.003 



87.708 - 90.000 
14.30 - 17.80 
0.00- 1.00 

0.1000-0.1125 
0.741 - 0.922 



flux value, before normalising the light curve (the uncertainty 
is accounted for separately, see Section [3). The full normalised 
light curve is shown in Figure |2l We evaluate the actual noise 
level per 512 s by measuring the dispersion about 1 h (7 expo- 
sures) bins and scaling it by V7, giving 8.9 x 10"^, compared 
to a photon noise level of 5.6 x 10 Possible factors contribut- 
ing to the difference include residual instrumental effects and the 
intrinsic variability of the star. 



3. Transit analysis 

A preliminary ephemeris (orbital period P and epoch Tq) was 
obtained by least-squares fitting of periodic, trapezoidal transits 
to the light curve after filtering out the out-of-transit variations 
using a 1-day baseline iterative non-linear filter. A more care- 
ful removal of the variability was then carried out by fitting a 
straight line to a light curve section lasting a little over one tran- 
sit duration before and after each transit. We also experimented 
with higher order polynomials, but they did not improve the dis- 
persion of the residuals, and do not change the results of the 
subsequent analysis. We folded the corrected segments of light 
curve using the preliminary period ephemeris, rebinned them in 
bins of 0.0003 in phase, and fitted the result to obtain preliminary 
estimates of the system scale a/Rs, the radius ratio Rp/Rs, the in- 
clination / and the linear limb-darkening coefficient m, where /?s 
is the star radius, Rp the planet radius and a the semi-major axis. 
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Fig. 3. Folded, binned light curve with the best-fit transit model. 



We opted to fit for u rather than fix it because reliable theoreti- 
cal limb-darkening coefficients are not currently available for the 
CoRoT bandpass. The ephemeris was then refined by fitting for 
the time of transit centre Tc for each individual transit event (fix- 
ing all other parameters) and fitting a linear relation to the Tc's. 
Finally, the light curve was folded again at the refined ephemeris 
and rebinned to perform a final fit for a/R^, Rp/Rs, i and u. 

At each stage, we use the formalism of iMandel & Agoll 
(I2OO2I) with quadratic limb darkening to generate model transit 
light curves and the Idl implementation Mpfit of the Levenberg- 
Marquart fitting algorithm, kindly provided by C. Markwart, to 
perform the fit. The period was fixed at the ephemeris value, 
and the epoch was also fixed except when fitting individual tran- 
sits. The eccentricity was assumed to be zero (the best fit to the 
radial velocity data is a circular orbit with an eccentricity uncer- 
tainty of 0.1, see Paper V). We also tried fitting the transits with 
a quadratic limb-darkening prescription, but this did not improve 
the fit, and therefore we reverted to linear limb-darkening. 

To evaluate the noise-induced uncertainties on the transit pa- 
rameters, including the effect of red noise, we used a 'correlated 
bootstrap' approach. The residuals from the global best fit were 
divided into bins lasting 1.12h (2/3 of the transit duration, or 
1/4 of the duration of the light curve segments used to cali- 
brate the out-of-transit variations around each transit), randomly 
shuffled, and added back to the fit before fitting the individual 
transits. Over-sampled and non over-sampled bins were shuffled 
separately. Each bin is shuffled whole, preserving the detailed 
time sampling of individual bins, so the procedure does not ac- 
count for the effect of small data gaps, but it does account for 
correlated noise on hour timescales, including the effect of star 
spots crossed by the planet. We used 100 realisations when fit- 
ting individual transits and 1000 when fitting the folded light 
curve. At each realisation, we also added a constant drawn from 
a Gaussian distribution with zero mean and standard deviation 
0.001 to the data, to account for the uncertainty in the contam- 
ination fraction. As the frequency distributions for each param- 
eter can be strongly non-Gaussian (see in Fig. |4]in the online 
material), we measure uncertainties as the interval away from 
the best-fit value where the frequency drops below e"'^^ times 
the maximum (if the distributions were Gaussian, this would be 
equivalent to the standard deviation). The results are reported 
in Table [T] For comparison, we also show in Fig. |4] the results 
of a standard bootstrap (which simply consists in swapping data 
points, and accounts for white noise only). Except for the limb- 



darkening coefficient, the two processes give similar results, in- 
dicating that red noise affects the other parameters' only slightly. 

To gain an insight into the effect of parameter-to-parameter 
correlations, we used a Bayesian approach. Reduced values 
were computed over a grid in 4-dimensional parameter space (/, 
a/Rs, Rp/Rs and u) about the best fit. A uniform grid in cos(0 
was used, which is equivalent to assuming an isotropic distribu- 
tion of inclinations. The grid was uniform in the other param- 
eters, i.e. no a priori information on these parameters was as- 
sumed (though u was restricted to the physical range 0-1). The 
X^'s were then converted to relative probabilities for each indi- 
vidual model using p cc exp(-x^ /2) and normalised. One can 
then obtain probability distributions for each parameter (shown 
on Fig.|4]i by marginalising over successive parameters. For each 
parameter, we report in Table[T]the interval over which the prob- 
ability is higher than e^'^^ times the maximum. This interval 
should be interpreted with some care: it is not a confidence in- 
terval in the frequentist sense, but rather an interval containing 
~ 68% of the posterior probability integrated over all other pa- 
rameter^ In the present case, this interval is very wide because 
the global minimum in the multi-dimensional surface is very 
narrow, but it is located at one end of a valley which widens 
significantly away from the minimum, as illustrated in Fig.|4] 

In Paper V, the transit parameters are combined with ground- 
based follow-up observations to give the stellar and plane- 
tary parameters, which we reproduce here for completeness: 
Teff = 6190 + 60 K, \Qgg = 4.41 + 0.05, M, = 1.16!°;°^M0, 
R^ = \.\1^_I^^Rq, age l+i;°GyrjMp = 0.72 + 0.08 Mjup and 
/?p = 1.19+JJ",5'5/?jup. Note that lClareti ( |2004 compute theoreti- 
cal linear limb-darkening coefficients in the range 0.55-0.6 for 
Teff = 6250 K and log^ - 4.5 in r', which is the standard band- 
pass closest to the peak of the CoRoT bandpass, though the later 
is much broader (300-lOOOnm). This is compatible with the 
value of 0.44 ± 0.15 obtained from the transit fit. 

We checked for transit timing variations in the O - C (ob- 
served minus computed) residuals from the refined ephemeris. 
The results are given in Figure |5] (online only). The third tran- 
sit in the time series shows a strong (> 500 s) deviation. To test 
whether this is a real timing variation, we repeated the individ- 
ual transit fits allowing ajR^ and Rp/Rs to vary as well as Tq, 
and found a clear correlations between the timing residuals and 
a/Rs (see bottom left panel of Figure |5]l, which points towards 
the effect of star spots or instrumental systematics rather than a 
real timing variation as the cause of the outlier. Closer inspection 
reveals that this transit contains small data gaps, and we interpret 
the deviation in measured timing and duration as an artefact of 
these gaps rather than a physical effect. We also checked for a 
secondary eclipse (removing the variability around phase 0.5 us- 
ing linear fits as done for the transits) but none is detected (as 
expected for this relatively low irradiation planet). 

4. Stellar rotation 

The light curve in Fig. |2] shows clear out-of-transit variability 
typical of a rotating, spotted photosphere. The semi-coherent na- 
ture of the light curve suggests that active regions on the star 
are evolving on timescales slightly longer than the rotation pe- 
riod. A given active region may therefore cross the visible hemi- 
sphere only 2 or 3 times in its lifetime. The form of the light 
curve evolves quickly, but partial coherency should persist for 
2 or 3 rotations, making auto-correlation an appropriate method 

' For an excellent discussion of Bayesian inferenc e and model s elec- 
tion and how it differs from frequentist methods, see lTrottj l l2008h . 
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Fig. 6. Weighted auto-correlation function of the out-of-transit 
Hght curve (gray Hne). We estimate the rotation period of the 
star using a Gaussian fit to the peak of the first side-lobe (dashed 
black line). 

to measure the stellar rotation period. We used an inverse vari- 
ance weighte d adaptation of th e Discr et e Correlation Function 
method of Edel son & KrolikI (1 19881) dCollier Cameron et al.1 
12008 ) to compute the Auto-Correlation Function, or ACF (see 
Fig.©. We fit a Gaussian to the first peak in the ACF to deduce 
a period of Pacf = 8.87 ± 1 . 12 d. The dominant source of uncer- 
tainty is the limited duration of the time-series. This is consistent 
with the V sin / of 6.4 + 1 .0 km/s reported in Paper V. 

5. Conclusions 

CoRoT-Exo-4b's orbital period is the secon d longest of any tran - 
siting planet to date, following HD 17 156b (iBarbieri et al.l2007l) . 
To our knowledge, it is only the second transiting system where 
the star's photometric rotation period ha s been measured pho - 
tometrically (the other being HD 189733, Henry & Winnll2008l) . 
There are other systems where the rotation perio d and planet 
orbital period are very sim ilar, particu larly t Boo dCatala et alj 
20071; iD onati et al.' '2008*) and XO-3 dJohns-Krull et al.l l2008t 



Winn et al. 2008) (the rotation periods are derived from Doppler 



imaging and v sin / respectively, rather than measured photomet- 
rically). In both cases, the period of the orbit and the stellar ro- 
tation period are not exactly the same, but is compatible with a 
synchronised outer envelope when allowing for differential sur- 
face rotation. Both planets, which orbit F-type stars like CoRoT- 
Exo-4, are massive and have very short periods (~ 3 d), and thus 
may induce strong enough tides in the star to have synchronised 
its convective envelope. 

A simple calculation, base d on the work of 
iDobbs-Dixon et al . (2004) and Jack son et all (12008). shows that 
CoRoT-Exo-4b, which has a smaller mass and longer period, 
would not exert a significant tidal torque on its host star unless 
one invokes an unphysi cally high rate of tid a l ener gy dissipation 
in the star However, iMcCuUough et al.l (l2008l) discuss the 
case of two other transiting systems containing F-type stars, 
XO-4 and HAT-P-6, both with orbital periods longer than 7 
days which are roughly half the (spectroscopically determined) 
stellar rotation period. They argue that this may indicate the 
existence of resonant interactions between the planet's orbit and 
its rotating host star, the absence of such resonances in systems 
containing cooler stars suggesting that Jupiter-mass planets 
can only interact effectively with stars with very shallow outer 
convective zones. 

Similarly, if some other factor brought the CoRoT-Exo-4 
system to the 1 : 1 ratio between stellar rotation and orbital pe- 
riod that is observed today, resonant interaction between star and 
planet may have maintained this ratio thereafter Any subsequent 
evolution of the star's rotation rate is likely to have been mod- 
est: 9 d is close to the peak of the rotation period distribution 
for F-type stars both during the TTauri phase and in the field 



(see e.g. Fig. 1 of lBarnesll2003h . The initial resonance may have 
occured naturally if the proto-planetary disk was truncated near 
the co-rotation radius, and the planets' migration halted close to 
the inner edge of the disk. This hypothesis is supported by the 
absence of detectable eccentricity (see Paper V). 

This system clearly warrants further observational and theo- 
retical investigation to pin down its tidal and rotational evolution 
status. For example, more detailed analysis of the out-of-transit 
light curve should enable the active region s on the stellar surf ace 
to be mapped in a time-resolved fashion ( Lanza et al.ll2007l) to 
search for signs of star-planet magnetic interaction. 
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Fig. 1. Correction of the hot pixel event in the blue channel. The 
light curve before correction is shown in grey, and after correc- 
tion in black. The vertical solid lines show the time limits of 
the correction section, with the exponential decay fit shown as 
the smooth red line. A small subsection, between vertical dashed 
lines, falls in-transit, and was not used to estimate the correction, 
though the correction was applied to it. The horizontal dashed 
line shows the reference flux level adopted before the hot pixel 
event. 
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Fig. 5. Transit timing residuals versus time (top panel) and com- 
pared to variations in transit duration (measured through the 
scale parameter ajRs) and depth (measured through the radius 
ratio Rp/R^. The apparent timing deviation of the third transit is 
accompanied by a change in duration, but not in depth. Closer in- 
spection of the light curve shows that this transit is more strongly 
aff'ected by data gaps than the others. 
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Fig. 4. Top: probability and frequency distributions obtained from the Bayesian and bootstrap methods respectively for the main 
parameters, ajR^, cos(i), u and Rp/Rs of the transit fit. The Bayesian method results are shown in grey, the correlated bootstrap in red 
and the non-correlated bootstrap in blue. Bottom: Two-dimensional probability distributions obtained from the Bayesian method for 
selected pairs of parameters, highlighting the correlations between parameters (particularly a/R^ and /). In each panel, the location 
of the best-fit model is marked by a red vertical arrow or cross. 



